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Abstract 



. Inspired by the bridge pioneered by Guerra among statistical mechanics and analytical 

mechanics on 1 + 1 continuous Euclidean space-time, we built a self-consistent method to 
r~| . solve for the thermodynamics of mean-field models, whose order parameters self average. 

We show the whole procedure by analyzing in full details the simplest test case, namely 
the Curie- Weiss model. Further we report some appHcations also to models whose order 
parameters do not self-average, by using the Sherrington-Kirkpatrick spin glass as a guide. 



> ' Introduction 
oo 

Mean field statistical mechanics of discrete systems is experiencing a massive increasing of in- 

^ ' terest for several reasons. Born as an infinite dimensional limit of a theoretical background for 

—I • ° 

condensed matter physics, mean field statistical mechanics become immediately appealing for its 
O ■ possibility of being solved (even though this happens exactly for really a few models [l4|), still 

j> I retaining several features of more realistic systems with finite dimensionality. 

Furthermore, and maybe nowadays, foremost, its range of applicability is continuously spread- 
ing such that, so far, it is one of the key tools for the investigation of several models far away 
5^ I from physics like biological or social networks (see for instance, respectively, [12] [13] and [5] |15|): 

all systems where the mean field nature of the description is not a limitation and whose rigorous 
or heuristic analysis was, in past decades, unimaginable. 

As a consequence the need for methods in statistical mechanics is one of the fundamental 
enquiries raised to theoretical physicists and mathematicians involved in the field. 

In this paper, inspired by a pioneering work of Francesco Guerra [8], we develop an alternative 
approach to standard statistical mechanics to solve for the thermodynamics of systems whose 
order parameters self-average. 

With the aim of presenting the theory also to readers who may not be experts in statistical 
mechanics, we apply our scheme to the simplest and most well known Curie- Weiss (CW) model, 
which we solve in full detail, for the sake of simplicity, linking our procedures with general 
statistical mechanics models via frequent remarks spread throughout the whole paper. 

As the largest interest is payed to complex systems, after the CW, we analyze the Sherrington- 
Kirkpatrick (SK) model, in the replica symmetric regime, subjected to an external field. 
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Guerra Action for mean field spin models 



Even though we will be interested in observable's behavior once the thermodynamic limit is 
taken, let us consider a large set of Ising spins cjj = ±1, i G (1, ...,N). Let us deal with a 
generic mean-field spin model, described by the Hamiltonian 

JV 

Hn{o-) = - ^ XijO-iO-j, (1) 

where Xij is a two body interaction matrix. The main quantity of interest in statistical mechanics 
is the infinite volume limit of the free energy = lim^^oo fN{P) = ^^^N^oo —P~^An{P), 
where An ((3) is the pressure and is related to the Hamiltonian via 



N 

We stress here (even though we will not deal with disordered systems in the first part of the 
work) that for the SK model it is usually expected to consider the quenched average of the free 
energy [8], however, without explicit expectation over the random coupling we mean its value 
X-almost surely in the sense of the first Borel-Cantelli lemma. 

It is useful to consider the one body interaction, of the same nature of Hamiltonian, that we call 
cavity field 

N 
i 

We define further a two parameters Boltzmannfaktor B{x, t) and a relative Gibbs measure {■)(x,t) 



as: 



BN{x,t) = e^p{9{t)HN + e{x)H'N) , (2) 
_ E.fi<r)iB{x,t)) 

where is a increasing function, vanishing at the origin, strictly dependent by the form of 
interaction. Eventually a magnetic field can be added in ([T|), and therefore in (I2|3p . 

We define the Guerra action (p{x,t) for a mean field model as the solution of the Hamilton- 
Jacobi differential equation 

dt'fN{x,t) + ^ {dx'fN{x,t)) + VN{x,t) = 0, (4) 

with suitable boundary condition. 

Furthermore the function t) = dx(p{x,t) satisfies 

dtUN{x, t) + un{x, t)dxUN{x, t) + dxVN{x^ t) = 0. (5) 

The Guerra action ipN{x,t) is related to the pressure of the model ^Ar(cj), in a way that will be 
specified later, case by case. 

Consequently even the potential function VM{x,t) expresses thermodynamical quantities of 
the case study (i.e. in CW and SK models we investigate, it turns out to be the self- averaging 
of the order parameters). 
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We will be interested throughout the paper in the region where V{x, t) = 0, when we can 
always solve our equations Q- In fact these problems are largely studied in literature far away 
from statistical mechanics [4]. In particular some Theorems, due to Lax [10], are helpful, since 
under certain hypothesis (that in a nutshell are related to the uniform convexity of the quantity 
\{dxip{x,t))'^), give the form of the unique solution of ([!]) and ([5]) (related by a derivation). 
Following Lax we can state the next 



Theorem 1. For a general differential problem 

dtip{x,t) + ^{dx(p{x,t)f 
if{x, 0) = h{x) 



m 



and 



dtu{x, t) + u{x, t)dxu{x, t) 
u(x,0) = g{x) 



m i 
on 



X (0,+cx)) 
x{t = 0}, 

i X (0, +oo) 
lx{t = 0}, 



(6) 



(7) 



where h{x) is Lipschitz-continous, and g{x) 
function y{x, t) : R x R+ M such that 



ip{x, t) 



h'(x) G C°° , it does exist and it is unique the 




+ Ky) 



+ h{y{x,t)) 



(8) 



y{x,t) 



(9) 



is the unique weak solution of and 

u{x, t) = 

L. 

is the unique weak solution of Furthermore, the function x y{x,t) is not- decreasing. 

It is worthwhile to remark that the choice of looking for weak solution (that arises naturally 
in the Lax's theorems) may look as redundant in our case, since we deal with physical quantities 
(in general smooth functions). Actually it prevents us from the eventual discontinuities of the 
solutions of jH]) and However, a strong solution is a weak solution too and there is no need 
to change the essence of the Theorem. 

Let us start applying this framework to the CW model. 



Mean field ferromagnet as a 1-dimensional fiuid 

The mean field ferromagnetic model is defined by the Hamiltonian 

^ N N 

^We stress however that the formalism we develop can still be applied to general constrained problems 
{V{x,t) 7^ 0), even though their resolutions can be prohibitive 
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It is easily seen that we have resumed in the Hamiltonian both the two body and one body 
interactioiil. Thus, choosing 6{a) = a, we can write the (rc, t)-dependent Boltzmannfaktor as 

I t ^ TV 

Bn{x, = exp — ^ (JiCTj + X ^ cJi 

Remark 1. When dealing with the ferromagnetic Boltzmannfaktor B^ix^t) above, classical sta- 
tistical mechanics is recovered of course, in the free field case, by setting t = (3 and x = 0. 

In the same way the averages {■)(^x,t) denoted by (.) whenever evaluated in the sense of 

statistical mechanics. 

A fundamental role is played by the magnetization m which we introduce as 

m = lim rriM = lim y^(Tj, (m) = lim ^ — £ — r-rr^- 

Let us denote UN{x,t) the 1-dimesional velocity field and if]\f{x,t) its dynamic potential (such 
that dx(p{x,t) = u{x,t)). Here the label N remembers us that the analogy is made with the 
CW model with finite size N (of course we are interested about the thermodynamic limit of the 
model). 

The Guerra action can be written as 

^N{x,t) = -^log exp (^^Nm% + xNiTiN^ = -^^(x, t) + O (^^^ , (10) 

Wn} 

i.e., the mean field CW pressure (up a minus sign) [3], where t stands for the inverse temperature 
P and X takes into account the external magnetic field h. 
Deriving (fTO]l we get 

UNix, t) = - {uIn) (x, t), (11) 

the mean value of the magnetization. So our analogy is now completed, and we can write a 
fluid equation as a transport equation for ui\i{x,t), plus an Hamilton-Jacobi (HJ) equation for 
(pN{x,t) and a continuity equation, defining the (purely fictitious) density function p{x,t). 

We notice that the Guerra action ip]y{x,t) satisfies an HJ equation where the potential 
function is the self- averaging of the magnetization. Indeed, since we have 

dt(pN{x,t) = {m%) , 

and 

dl2ipN{x,t) = dxUN{x,t) = -dx {ruN) {x,t) = -N{{m%) {x,t) - {rriNf ix,t)), 

we can easily choose the external pressure for the fluid, that appears as a potential in the HJ 
equation, as 

Vn{x, t) = ^[ {ml) (x, t) - {mN? (x, t)) , (12) 



^In ferromagnet the cavity field coincides with the external field 
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and we have also 

-^dlipN{x,t) = VN{x,t). (13) 

Finally, computing 

(^Ar(3;,0) = —An{x,Q) = — log 2 — log cosh X, (14) 
we can build the differential problem for our hydrodynamical potential (piy{x,t): 



dtipN{x,t) + ^{dxipN{x,t)f + VN{x,t) = in M X (0,+oo) 
(Pn{x,0) = — log 2 — log cosh X on M x {t = 0}. 



(15) 



Remark 2. We stress that by choosing as a boundary a general point on x but t = (as we did 
in eg. [14])), we implicitly skipped the evaluation of the two body interaction which is, usually, the 
hard core of the statistical mechanics calculations as the one body problem trivially factorizes. 



Eq. (fT5|l is just the Hamilton-Jacobi equation for a mechanical 1-dimensional system, with 
time-dependent interactions. We can write it in a more suggestive way, for exalting our hydro- 
dynamical analogy. Indeed, bearing in mind (|13p . we have 

r dtfN{x,t) + ^{dxipNix,t)f-jj^dl(pNix,t) = inRx(0,+oo) 

I v9Ar(x, 0) = — log 2 - log coshx on R X = 0}. ^ ' 

This equation is more interesting than the first one, for several reasons. At first it is closed 
with respect to the unknown functioiil. Furthermore it has a clear physical and mathematical 
meaning: Indeed the presence of a dissipative term suggests the typical viscous fluid behavior, 
where friction acts against the motion. The smallness of this term (that appears with a factor 
N-^) acts mollifier for our differential problem. It may appear even clearer by investigating 
the equation for UAr(x,t). Deriving with respect to x eq. ffTB]) (and using standard results about 
for the order of derivation) we obtain 

dtUN{x,t) + UN{x,t)d^UNix,t) - ^d'^UN{x,t) = inMx(0, +00) , . 

UN {x, 0) = -tanh{x) on M x {t = 0}. 

This is a viscous Burgers' equation, i.e. a very simple Navier-Stokes equation in one dimension. 
Here the mollifier term is more incisive, since, as we will see soon, when it vanishes (i.e. in 
thermodynamic limit), it induces the spontaneous Z2 symmetry breaking of statistical mechanics 
by making the solution u{x, t) (i.e. the magnetization) not regular in the whole (x, t) half-plane. 

Lastly let us derive the continuity equation that should complete our formal hydrodynamical 
analogy for the ferromagnetic model. We stress that it does not carry any further information 
about the model, as it is all contained in (fT6]l and (fTTll ). From the continuity equation we get 

dtPN{x,t) +UN{x,t)dxPN{x,t) = -pN{x,t)dxUNix,t) 

= pN{x,t)2NVN{x,t). 

Writing 

DN{x,t) = dt +UN{x,t)da: = —, (18) 

as 



''This is actually a feature of the ferromagnets. For instance it is easily seen that it is not trivially closed for 
the SK pressure because every derivation involves different overlap combination [2]. 
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the differential operator along the stream lines, we obtain the equation for 

Dn{x, t)pN{x, t) = 2NVn{x, t)pN{x, t), (19) 



solved by 
that is 



=pjv(0,0)e2^^'^"^(^(^)'*(^)) (20) 

Pn{x, = ^ exp [Ntm% + Nxtun] =^Ar(2t,x). (21) 

Resuming, mean field ferromagnets of finite size is completely equivalent to the 1-dimensional 
viscous fluid described by equations 

dtUN{x, t) + Un{x, t)dxUN{x, t) - 27^5^2 WAT (x, t) = 

DN{x,t)pN{x,t) = 2NVN{x,t)pN{x,t), 

and in thermodynamic limit, to an Eulerian fluid, such that 

J dtu{x,t) + u{x,t)dxu{x,t) = 
\ p{x,t)-^D{x,t)p{x,t) = 0. 

We would like now to link the finite dimensional model with its thermodynamic limit, i.e. the 
viscous fluid with the inviscid one. It is consequently useful to study the free problem 

dtifix, t) + ^{dM^, t)f =0 in M X (0, +oo) 
V9(x,0) = — log2 — logcoshx on M X {t = 0}, 

and 

dtu{x,t) + u{x,t)dxu{x,t) = in M X (0, +CX)) , , 

u(x,0) = -tanhx on M x = 0}. ' ^ 

With this purpose we can use Theorem [D 

Remark 3. We stress that via Theorem (OP changing the boundary condition is equivalent to 
modify the nature of the spin variables in the ferromagnetic model. Since the condition on H is 
Lip schitz- continuity, such a theorem is valid for every distribution of spin variables with compact 
support, but not for example for Gaussian ones (at least trivially). We let for future works further 
investigations JB]. Hereafter anyway we will deal with only dichotomic variables. 

With h{y) = — log 2 — log cosh y, y = x — tu{x,t) (given by jH])), we find 

if{x, t) = -u{x, i)^ — log 2 — log cosh (x — tu{x, t)) , 

and bearing in mind if = —A and u = — (m), by setting t = /3 and x = h, we gain the usual free 
energy for mean field ferromagnet 

111 ( B (m)"^ 1 
/(/3, h) = --A{(5, h) = -ifiP, h) = -l - log cosh /3 (/i + (m)) - log 2 I , 

where of course (m) is the limiting value for the magnetization, as we are going to show. We only 
have to prove convergence for differential problems (fT6]l and (fTTl) to the free ones, respectively 
([22]) and (f23ll . Let us start with the former by stating the following 
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Theorem 2. The function 



fNix,t) = -—log 



+00 



dy 



-00 V 27r 

solves the differential problem (E^, and it is 



exp —N 



{x - yf 
2t 



log 2 — log cosh y 



(24) 



\(pN{x,t) - f{x,t)\ < 0{^). 

Proof In order to find a solution of ([22]l . we puiEl 

After a few calculations 

dt(j)N{x,t) -- 



(25) 



]-N(j)N{x,t){dx(pN{x,t))'^ - l;(t>N{x,t)dl2(pN{x,t) 



^dl2^N{x,t), 



(26) 



we see that (l){x, t) solves the heat equation with conductivity (and a suitable boundary 
condition) : 



dt(j)Nix,t) - ^dl24>Nix,t) = in M X (0,+oo) 
(t>N{x,0) = 2-^cosh-^2; 



on M X {t = 0}. 



(27) 



The unique bounded solution of (|27|l is 



4>N{x,t) 



+00 



dy 



exp -N 



1-00 \/27r 

and, bearing in mind 997V = — log i;^>Ar, we have 



(x - yf 
2t 



log 2 — log cosh y 



(pN{x,t) 



N 



log 



+00 



dy 



exp 



-N 



yf 



2t 



log 2 — log cosh y 



We notice that, since the uniqueness of the minimum of the function in the exponent (allowed 
by Theorem ([I|)), we easily get fN^f when — > 00. 

Finally bounds on the error can be made via standard techniques. □ 

We must now prove an analogue result for the velocity field u{x,t). Since the equations 
for ip{x,t) and u{x,t) are trivially related by a derivation, it is clear that UN{x,t) — > u{x,t) 
uniformly in the thermodynamic limit. Anyway for the sake of completeness (and as a guide for 
testing other models) we state the following 



*This is usually known as the Cole-Hopf transform i4j. 
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Theorem 3. The function 



/-» ^ ( - - log 2 - log ™h y] 



solves the differential problem Ii23\} and it is 



\uNix,t)-u{x,t)\<0(—^]. (29) 



Proof The (1281) is easily obtained by direct derivation of ifN in ([SI 
Again the bound on the error is made via standard techniques. □ 
Finally we can state the subsequent 

Corollary 1. It is VN{x,t) < O(^) a. e.. 

Proof For the two previous theorems we have 

ipN(.X,t) = (p{x,t) + O(^), 

thus ^ 

dtifN = dtip + 0{ — ) 



and 



and therefore, using the Hamilon-Jacobi equation (fTSl) for ^pN, we find 

that implies the thesis. □ 

What we meant for "a.e." is actually the whole (x,t) positive half-plane, but the line defined 
by (x = 0, i > 1) as will be well explained in the next section. 

Shock waves and spontaneous symmetry breaking 

In this section we study more deeply the properties of equation ((231) . This is an inviscid Burgers' 
equation, and again we can have a representation of solutions as characteristics j^. We get 

u(x, = — tanh(x — ii(x, (30) 

i.e the well known self consistence relation for the CW model, with trajectories (parameterized 
by s G M) 

■ * = ' (31) 



X = Xq — s tanhxQ. 
We can immediately state the subsequent 
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Proposition 1. In the region of the plane defined by 



> -y^t{t-l) + arctanh ( J- ] for xq > 



and 



< -^/t(P^ + arctanh { J for xq < 



t- 1 



trajectories I131\) have no intersection points. 

Remark 4. This last statement defines the onset of ergodicity breaking in the statistical mechan- 
ics of the CW model. 

Proof Set for instance xq > 0. 

Once fixed s = s let us investigate the position at time s as a function of the starting point 
Xq. We have 

x{xo) = xq — stanha^o- 

If x{xo) is monotone with respect to xq, then Vxq G K 3!x(t), i.e for every starting point there 
is an unique position at time t. In other words, two trajectories born in different points of the 
boundary cannot, at the same time, assume the same position (do not intersect). Hence we have 

x'{xo) = 1 — s(l — tanh^ xq) > Vxq, 

only if 

1 

s < 



1 — tanh^ Xq ' 



as 1 — tanh^ Xq always belongs to [0, 1]. The last formula implies 

Xq > arctanh 



t - 1 



t 

and bearing in mind the form of trajectories ([3T]) we get 



t-l 



X > arctanh y^J j - \/t{t - 1). (32) 

The proof is analogue for xq < 0. □ 

We must notice that the previous proposition gives the region of the (x, t) plane in which the 
invertibility of the motion fails. On the other hand, every trajectory has its end point at the 
intersection with the x-axes, or are all merged in a unique line, that is (x = 0, t > 1). 

More rigorously, the curve (x = 0, t > 1) is a discontinuity line for our solution, since it is 
easily seen that every point of such a line is an intersection point of the trajectories ([3T]l . Also 
we can get by ((30]) with direct calculation 

1-v? 

9Xx,t) = - ^^^^^_^,^ <0, (33) 
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i.e. the velocity field is strictly decreasing along x directioiil. 

Now we name m+ the limiting value from positive rc, and m_ the one from negative x, and state 
the following 

Proposition 2. It is < u_ = —u^ < for a.e. t > 1. 

Proof The curve of discontinuity can be parameterized as 

t > 1 
X = 0, 

so has zero speed. We have that V t > does exist a neighbors / of (x = 0) such that u{x,t) is 
smooth on I. Thus, since we know that our u{x, t) is an integral solution, we can use Rankine- 
Huginiot condition [4] to state 

2 2 
= U_. 

Since for ([331) it has to be uj^ < u_ the assert is proven. □ 

Remark 5. We stress that the relation = u^, in this context, mirrors the spin-flip symmetry 
shared by the two minima of the CW model in the broken ergodicity phase, i.e. \ + {m)\ = \ — {m)\. 

It follows that (x = 0,t > 1) is a shock line for the Burgers' equation (f23]l . 
On the other hand, of course, x = is an equilibrium point for the system, since we have 
that both position and velocity are zero. 

Remark 6. This property is translated in statistical mechanics to the trivial case of CW model 
without neither a vanishing external field, such that spontaneous magnetization can never happen. 

We can use it for exploring the well known mechanism of spontaneous symmetry breaking. 
With this purpose, let us move on a family of straight lines of equation 



X = et — e. 



We have infinitely many lines, all converging in (0,1), that intersect the x-axes in — e. Let us 
choose for example e > 0, and perform the limit of u{x, t) on the shock line taking the value of 
u{x,t) by these, and then sending e — > 0. Since — e is negative, the intersection point with t = 
is approaching from the left (x~), meanwhile the limit of u is taken from the right In 
the same way we have that when the intersection point approach to zero from right (x"*"), the 
limiting value of u is taken from left (u-). 

Remark 7. In our analogy with statistical mechanics one can make the substitution n(x, t) = 
— (m) {h,l3), and t = (3, x = h(3, getting the spontaneous symmetry breaking mechanism, in such 
a way that limf^^Q±{m){h, P) = . 

^This is a particular case of a more general property of the Lax-Oleink solution |10| . named entropy condition^ 
that ensures u{x, t) never increases along x. We won't give the general form, that is redundant in this contest, 
but can be very useful in studying generalized ferromagnet [6]. 
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Conservation laws 

We can rewrite the (fT5]) from a mechanical point of view as 

dtifN{x,t) + HN{dxipNix,t),X,t) 

and the Hamiltonian function reads off as d 



HNid^ifNix, t),x, t) = + Vn{x, t). (34) 

Hamilton equations are nothing but characteristics of equation (fT5]l : 

X = Uiy{x,t) 

' = ^ (35) 

p = -UN{x,t)dxUN{x,t) - dxVN{x,t) 

E = -UN{x,t)dx{dtipN{x,t)) - dtVN{x,t), 

the latter two equations express the conservation laws for momentum and energy for our system, 
and can be written in form of streaming equations as 

DNUN{x,t) = -d^VNix,t) 
DN{dt(pN{x,t) = -dtVN{x,t). 

Since in thermodynamic limit the system approaches a free one, bearing in mind that U]y{x, t) = 
— {uin) and dt(pN{x,t) = — ^ {rnj^'j, so Dn = dt — {mj\[) dx, for ^ cxd we conclude 



Dn {rriN) = 
DN{m%) = 0, 



(36) 



I.e. 



{m%)-3{mN){m%) + 2{mMf = O(^) 
(m2)^)-2(m)(m3> + 2(m)2(m2) = O(^). ^ ^ 

We have from Corollary [1] that (^m?'j = (m)^ + 0{j^) everywhere but on the line {x = 0,t > 
1), where anyway (m) = 0. It is possible to write down a relation that follows from energy 
conservation: where the potential vanishes, using momentum conservation, giving {|m'^ 
0{jf), we get 

(m^> - {m'f = 0(1). 

Otherwise when the potential is different from zerc[l| we have (m) = 0, thus the previous formula 
is still valid, and it holds in all the (x, t) half-plane. 

Remark 8. This is of course a Ghirlanda-Guerra relation l?^ for the CW model (i.e. it expresses 
self-averaging of the internal energy density). As a counterpart, the bare momentum conservation 
implies the first Aizenman-Contucci JJ] relation for (m^). 

Remark 9. It is interesting to remark that the orbits of the Mother groups of the theory coincide 
with the streaming lines of our fluid, and conservation laws along these lines give well known 
identities in the statistical mechanics of the model. 



here we namep our velocity u, i.e. the velocity field coincides with the generalized time dependent momentum 
'^Anyway it is a zero measure set. 
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The Replica Symmetric phase of the Sherringon-Kirkpatrick model 



Despite the main goal when dealing with complex systems is a clear scenario of the Broken 
Replica Phase, which, in our languages translates into solving viscous problems as VAr(x,i) / 
(and it is posted to future investigations), a detailed analysis of the replica symmetric regime is 
however immediate within this framework, as pioneered in [8]. 



The Sherrington-Kirkpatrick Hamiltonian is given by 



where Jij are i.i.d centered Gaussian variables, with E[Jjj] = and IE[Jj^] = 1. 
Following [8] we introduce the partition function 

Zn{x, i) = ^ exp I ^ JijCTiCTj + Vx'Y^ JiCTi + ^ fTj 

Accordingly with the normalization factor \/\fN of the model, we choose 6{a) = y^; it is 
important to stress that differently to the ferromagnetic model, the cavity field with strength y/x 
does not coincide with the magnetic field /i, that entries in the Boltzmannfaktor as an external 
parameter. Thus our results will hold for every value of h. 

The main difference, when introducing thermodynamical quantities (as the free energy) is in 
an overall average over the random quenched couplings encoded in the interaction matrix. In 
this sense the averages (.) now stand both for the Boltzmann averages (denoted by uj hereafter 
when dealing with a single set of phase space configuration, = ujxujx ...>iuj when dealing with 
several replicas of the system) and for the averages over the coupling (denoted by E hereafter), 
such that (.) = E0(.). 

The Guerra action for the SK model reads off as 

ipN{x,t) = 2An -l^- X. (38) 

So it has, once introduced the two replica overlap as qi2 = N^^ af''* af'\ 

dtipN = 2dtAN-\ = -\{ql2) (39) 
d^^N = 2d^AM-l = -(gi2>. (40) 

Mirroring the mean field ferromagnet, also in this glass model the interaction factorizes at f = 0, 
and, once set E„ = —k= f ^ dge~^ , we have 

ipN{x, 0) = 2Af/^{x, 0) - X = 2 log 2 + 2Eg log cosh(^/i + g^/^) - x. 
The last formula, together with ([39l [40]) allows to build the HJ equation for (pN{x,t) 



dtipN{x,t) + i(5^(/7Ar(x,t))2 + VN{x,t) =0 in M X (0,+oo) 

ipNix, 0) = 2 log 2 + 2Eg log cosh(/3/i + g^/^) -x on R x {t = 0}, 



(41) 
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with 

VN{x,t) = ^[{ql,)-{quf). (42) 

In complete generality, this is an equation more complicated than the ferromagnetic one: 
Reflecting the complex structure of the RSB phase, the closure of the equation can be obtained 
only via cumulant expansions of the overlaps in terms of higher order correlation functions [2] , 
i.e the potential has no trivial expression in terms of (p]\f derivatives. We will study this equation 
in the Replica Symmetric phase, that is where, in the (x, t, h) domain, lim^ = 0. 

The velocity field, accordingly with ([40]) . is 

UN{x,t) = - (qu) {x,t) 

and satisfies the transport equation 



dtUN{x, t) + un{x, t)dxUN{x, t) + dxVN{x, t) = in M X (0, +oo) 
un{x, 0) = -Eg tanh2(/3/i + g^) on M x {t = 0}. 



(43) 



Remark 10. We stress that naturally in our approach the hyperbolic tangent of the CW model 
has been mapped into the squared hyperbolic tangent in the SK case, exactly as it happens in 
statistical mechanics, reflecting the role of the overlap as a proper order parameter with respect 
to the magnetization. 

Replica symmetry apart, the characteristic trajectories of ([43]l are not in general straight 
lines, because of the presence of the potential. We can give an expression for them: 



t = s 



X = xq — sEgtanh {j3h + q^^xq) — ds' dxVN{x{s),t{s)). 



(44) 



and solving for u 



U]\[{x,t) = — Eg tanh^(/?/i + g(A/xo(x, t)) — / dsdxVN{x{s), s), (45) 



where we get xo{x,t) inverting the second among (|44l) . 

This is the analogous of the Guerra sum rule for the order parameter stating that the 
difference among the true order parameter and the RS one is the line integral of the x derivative 
of V/v along trajectories. 



Reducing our attention to the RS phase of the model, we get the free HJ equation 

dt(pRs{x,t) + ^{dx(pRs{x,t)f = in M X (0,+oo) 

ipjisix, 0) = 2 log 2 + 2Eg log cosh(/3/i + g^/^) -x on M x {t = 0}, 



and Burger's equation 



dtUj^{x, t) + un{x, t)dxU]\f{x, t) = in M X (0, +oo) 
UNix, 0) = -Eg tanh2(/3/i + g^) on M x {i = 0}. 



(46) 



(47) 



'Actually Guerra relation may be obtained thought an integration of l|45p. 
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We are now in perfect agreement with the hypothesis of Theorem ([T|). Therefore we can write 
the Rephca- Symmetric Guerra action, in the thermodynamic limit, as 



t f X — y{x, t) 



+ log 2 + Eg log cosh(/?/i + gyjy{x,t)) - y{x, t) 



(48) 



and, naming the velocity field of the free problem —q{x,t), we trivially get from (l45]l the self- 
consistence equation 

q{x,t) = Kg tanh^(/?/i + gy^ x + tq{x, t)), 



and the trajectories are 



t 

x 



xo — sEg tanh^(/3/i + g^yxo). 



(49) 



(50) 



Remark 11. We stress that eq. is exactly the self-consistent equation for the SK model 

order parameter in the replica symmetric ansatz. 

Furthermore the minimization point y(x, t) is usually given by 

y{x^ t) = X + tq{x, t). 
Proposition 3. For values oft, x and (5h such that 



tE„ 



cosh {(5h + g^Jx + qt) 



1 2 

< - + -tKn 

— Q Q y 



cosh {(5h + gyjx + qt) 



(51) 



trajectories ^5U\) have no intersection points. In particular the whole region with x > and t > 
is included in i51\). 



Remark 12. We notice that the h51\) gives the form of the caustics for the {x,t) motion, i.e. 



cosh {(3h + g^/aT+lji) 



1 2 
3 3 ^ 



cosh {Ph + g^/x + qt) 



and in this sense completes the theorem given in fSj. 

Proof The procedure is just the same used in Proposition [H Starting from ([50]l . we put 

x{xq) = xq — tKg tanh^ {f3h + gy/x^) , 

i.e. the position depending by initial data, and let's study its monotony. Given the trajectories, 
it is clear that, whereas there is no intersection, x{xq) must be increasing, thus 

dxoxixo) = 1 - tEgdxo tanh^ {ph + g^/x^) > 0, 

(of course we can swap derivatives and Gaussian integral). So we have 

tEgdxo tanh^ (/3/i + g./x^) < 1. (52) 
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Now 



Egdxo tanh^ {f3h + g^/xo) 



Xq 
1 



tanh + g^/xo) 
cosh^ + g-^/iro) 

tanh (/3/i + g^/xQ) 
cosh^ + g^/x(j) 



d, 



E„ 



3E„ 



cosh^ [[3h + flfy^) 



cosh + g^/xo) 



2E„ 



2E„ 



tanh^ iyfih + 
cosh^ (yj3h + g^/xo) 



cosh (/3/i + gy^) 



where we have used the well known formula for Gaussian expectation Eg [gF{g)\ = Eg \dgF{g)\. 
At this point, putting the last expression in ([52]l we gain the (f5T]l . □ 

We can finally give the form of the Sherrington-Kirkpatrick solution for the pressure of the 
model [9][Il]. It is 

ARsm = Ars{^,P^) 



g2 

log 2 + Eg log cosh(/3/i + (7/371) + ^ (1 - qf 



(53) 



Conservation laws 

In the same way we did for the CW model, we can get relation among overlap from momentum 
and energy conservation laws, holding in RS regime. It is remarkable that the vanishing, in 
thermodynamic limit, of an overlap polynomial is associated to a Nother streaming of mechanical 
quantities. 

With the aim of deepen this last paragraph, let us stating the following 
Lemma 1. Given F{a^...a^) as a smooth, well behaved function of s replicas, we have 



N / 



a<b 



?a6 ~ "5 



Qa,s+1 "I 7^ Qs+l,s+2 



The proof of this lemma works via a long and direct calculation, and we will not report it 
here 

We stress that the linearity of D implies all our relations approach zero as O {1/N). 

We have, in general, that conservation laws for momentum and energy are given by the 
streaming equation 



Dn {qi2) 



-5^VAr(x,t) 
-2dtVN{x,t). 



(54) 
(55) 
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Of course in RS phase, the right hand term of (f54ll and ([55]) vanishes when ^ oo. Although 
such an approach does not give any information about the way they vanish in thermodynamic 
limit (we can always write it is o(l)), we can write down our relation without problems, since 
the presence of D forces them to be at least 0{1/N). Explicitly we get 

N {qi2 - '^Quqh + 3912934) - ^ {Q12) {qi2 - 4gi2923 + 3912934) = 0(1) 

^(912 - 4912923 + 3912^34) - ^(^2) (912 - 4912923 + 3912934) = 

i.e. conservation of momentum and energy along the streaming lines of the systems (or along 
free trajectories (fSOll ) implies that in the RS regime 

(9?2 - 49l29i3 + 39l29i4)(^,t) - {<ll2)(^^t) {112 - 4912923 + 39l2934)(^ J) < 
(912 - 49?2923 + 39l2934)(^,j) - (9l2)(a;,t) {<ll2 " 4912923 + 39l2934)(^_j) < 

Combining the previous results, we get a third relation 

(9^2 - ^<ii2<i2Z + 39?29i4)(^,t) - {Qi2? (9?2 " 4912923 + 39i2934)(^_t) < O j , (56) 

which, in particular we find physically meaningful, when setting x = and t = P'^, because 
the replica symmetric assumption on the vanishing of the potential is clearly reflected into the 
overlap labels in the last identity. 

If now we neglect the magnetic field {h = 0), as we are in the replica symmetric regime, the 
gauge symmetry holds such that the SK Hamiltonian is left invariant under the transformation 
(T — > (7CT, CT being a dichotomic variable out from the A^-spin Boltzmann average. Matching [2] 
and [8] in fact it is straightforward to check that gauging the energy conservation we get (again 
we stress that it holds only at h = 0, and obviously at t = e x = 0) 

(1 - {cih)) {qt2 - 49?29i3 + 39?29i4) < o 

and consequently 

(9^2 - 49?2923 + 39?2934) - O 

obtaining the well known relation constraining overlaps |ljt2j- 

Conclusions and outlook 

In this work we built a self-consistent method to solve for the thermodynamics of mean field 
systems, encoded by self- averaging order parameters. 

Such a method minimally relies on statistical mechanics, essentially just on the boundary 
conditions of our partial differential equations, and however, involves just straightforward one- 
body problems. 

Within our approach, that we tested on the Curie- Weiss prototype, we obtained the explicit 
expression for the free energy as a solution of an Hamilton-Jacobi equation defined on a 1 -|- 1 
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Euclidean space time, whose velocity field obeys a suitably defined Burger's equation in the same 
space. 

The critical line defining ergodicity breaking is obtained as a shock wave for a properly defined 
Cauchy problem. The behavior of the magnetization, thought of as this velocity field, both in 
the ergodic and in the broken ergodicity phases have also been obtained rigorously. 

As instruments involved in our derivation, we obtained rigorously also the existence of the 
thermodynamic limit for the free energy and the self- averaging of the order parameter. 

Despite the problems in relating conserved quantities and discrete symmetries, in our con- 
tinuous framework, Noether theory is straightforwardly applicable and gives the well known 
factorization of the momenta of the order parameter, as expected, being the Curie- Weiss a mean 
field model. 

Furthermore we applied the method even to the replica symmetric phase of the Sherrington- 
Kirkpatrick model, founding full agreement with Guerra's results and stressing other points as 
the study of the caustics (which shares some similarities with the AT line) and the study of the 
symmetries, which turn out to be polynomial identities, typical of complex systems (first of all 
the Aizenman-Contucci relations). 

We emphasize that, actually, we believe the method working for a large range of models 
{i.e. with general interacting variables as spherical spins, etc), several interacting spins as P-spin 
models, etc.). However, of course, it is still not enlarged to cover the case of not self- averaging 
order parameters (which is mathematically challenge even with already structured methods). 

Furthermore a certain interest should be payed trying to apply the method to finite dimen- 
sional problems. 

We plan to report soon on these topics and their possible applications. 
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